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SUMMARY 


An analysis utilizing laminated plate theory is developed to 
calculate the strain energy release rate associated with local 
delaminations originating at off-axis, angle ply, matrix cracks in 
laminates subjected to uniaxial loads. The analysis includes the 
contribution of residual thermal and moisture stresses to the strain 
energy released. Examples are calculated for the strain energy 
release rate associated with local delaminations originating at 90 
degree and angle ply (non-90 degree) matrix ply cracks in glass 
epoxy and graphite epoxy laminates. The solution developed may be 
used to assess the relative contribution of mechanical, residual 
thermal, and moisture stresses on the strain energy release rate for 
local delamination for a variety of layups and materials. 


KEYWORDS: residual thermal stress, hygroscopic stress, moisture, 
strain energy release rate, delamination, matrix crack 
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NOMENCLATURE 


A Delaminated area 
a Delamination length 
C Compliance 

E Modulus of elasticity 

G Strain energy release rate 

h Ply thickness 

K Number of plies in a region 
I Length 

m Number of local delaminations 
N Number of plies in uncracked laminate 
n Number of cracked plies 

{N} Vector of force resultants 

_P Applied axial mechanical load 

[q] Transformed reduced stiffness matrix 

t Thickness of a region 

U Strain energy 

u Strain energy density 

V Volume 

W Work 

w Width 

ai Lamina coefficient of thermal expansion in 
the fiber direction 

ol 2 Lamina coefficient of thermal expansion 
transverse to the fiber direction 
{a}k Lamina coefficients of thermal expansion in 
the laminate coordinate system 
{a} Laminate coefficients of thermal expansion 
Pi Lamina coefficient of moisture expansion in 
the fiber direction 

P2 Lamina coefficient of moisture expansion 
transverse to the fiber direction 
(P}k Lamina coefficients of moisture expansion 
_ in the laminate coordinate system 
jpj Laminate coefficients of moisture expansion 

AH Percentage moisture weight gain 
AT Temperature differential from cure 
temperature to test temperature 


2 



eij 

{£} 

{e°} 

{k} 

Oij 

{o} 

"Hxy.x 

v X y 


Strain tensor 

Vector of in-plane strain components 
Vector of mid-plane strain components 
Vector of curvatures 
Stress tensor 

Vector of in-plane stress components 
Laminate coefficient of mutual influence of 
the second kind 
Laminate Poisson's ratio 


Subscripts and superscripts 


k Ply number 

LD Locally delaminated region 
LAM Laminated region 

M Mechanical Loading 

T Thermal Loading 
H Hygroscopic loading 
x,y,z Laminate coordinates 


INTRODUCTION 

One common source of delamination in composite laminates 
are matrix cracks that form in the off-axis plies oriented at a non- 
zero degree angle to the applied load. If many similar plies are 
grouped together in the laminate, matrix cracks may form under a 
monotonically increasing load and initiate local delaminations [1]. 
More commonly, matrix cracks cause delaminations to form under 
cyclic loads, even if the laminate stacking sequence contains angle 
plies with only single ply thicknesses [2,3]. Hence, in order to 
predict the damage tolerance and fatigue behavior of these 
laminates, solutions are needed for the critical parameters that 
govern the formation of delamination from matrix cracks. 

Previous studies have contrasted the difference between local 
delaminations induced by matrix cracks and the more commonly 
studied edge delaminations resulting from mismatch in lamina and 
laminate properties [4-8]. In several studies, the strain energy 
release rate associated with these localized delaminations has been 
shown to be a useful parameter in predicting the onset and 
accumulation of local delaminations through the thickness in flat, 
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tapered, and curved laminates [3,9,10]. In references 3&9, 
predictions were accomplished using a simple closed form equation 
for G that was derived using laminated plate theory. The purpose of 
this investigation was to develop a G analysis based on laminated 
plate theory that includes the contribution of the residual thermal 
and moisture stresses to the strain energy release rate associated 
with local delaminations from matrix cracks in laminates subjected 
to uniaxial loading. 

STRAIN ENERGY RELEASE RATE ANALYSIS 


Background 

Figure 1 shows a simplified model of a composite laminate 
containing delaminations growing from a matrix crack. The 
composite gage length, I, is divided into a locally delaminated region 
(LD) whose length is the delamination length, a, and a laminated 
region (LAM) of length l-a. In reference 4, an equation for the strain 
energy release rate associated with local delamination growth was 
derived using the classical definition for G as 


2 dA 


( 1 ) 


where C is the laminate compliance and A is the surface area 
created by the local delamination. A series solution was developed 
assuming the composite displacements are the sum of the 
displacements in the laminated and locally delaminated regions, and 
the load carried by the two regions is equal to the applied load, P. 
This series solution resulted in the following expression 


2mW 2 uldElD tLAM^LAM/ 

where P is the applied axial mechanical load on the laminate, w is 
the laminate width, and m is the number of delaminations growing 
from the matrix crack. This parameter, m, would have a value of 2 if 
the cracked off-axis ply is in the interior of the laminate, and would 
have a value of 1 if the cracked off-axis ply is a surface ply. 
Furthermore, tj_AM is the laminate thickness, Elam is the laminate 
modulus (as calculated from laminated plate theory [11]), tLD is the 
thickness of the locally delaminated region that carries the load 
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(i.e., the laminate thickness minus the thickness of the cracked off- 
axis plies), and Eld is the modulus of the locally delaminated region 
(as calculated from laminated plate theory). Although figure 1 
depicts a local delamination growing from a 90 degree matrix crack, 
eq (2) has also been shown to be valid for local delaminations 
extending from angle ply matrix cracks where the cracked ply is 
oriented at an arbitrary, non-90 degree, angle to the load axis [8]. In 
equation 2, G depends on the layup and the location of the 
delaminated interfaces, which is determined by the through- 
thickness location of the cracked angle ply in the laminate, but does 
not depend on delamination size. 

In reference 6, a 3D finite element analysis was utilized to 
determine the strain energy release rate for local delamination in 
laminates containing 90 degree matrix cracks. The strain energy 
release rate for local delamination increased initially with the 
delamination length from the matrix crack but eventually approached 
a constant level asymptotically at a delamination length of 3 to 4 
ply thickness from the matrix crack. For a given delamination length, 
the 3D analysis yielded plots of G across the laminate width. The 
asymptotic values of total G for local delamination calculated near 
the free edge agreed well with eq (2). The agreement was found near 
the free edge instead of in the interior of the laminated width 
because the modulus terms in eq (2) were calculated using laminated 
plate theory and not beam theory. Hence, as was shown in ref [6], eq 
(2) is sensitive to the mismatch in Poisson contraction of the 
original laminate and the locally delaminated region. 

Analysis Derivation 

Although eq (2) has been verified using the results of a 3D 
finite element analysis [6], it includes only the contribution of the 
mechanical load to the strain energy released as local delaminations 
grow. Residual thermal stresses and the stresses due to absorbed 
moisture may also contribute to the strain energy that is released 
as the delamination grows. Therefore, a new analysis was developed 
to account for these thermal and moisture influences. 

Beginning with the definition of strain energy release rate, 

G =dW . dU. (3) 

dA dA 

where W is the work performed by the applied mechanical load and U 
is the strain energy released as the crack extends, an expression 
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may be derived for G associated with local delaminations from 
matrix cracks that includes the influence of residual thermal and 
moisture stresses. The work performed on the laminate by the 
applied mechanical load is equal to the summation of the work 
performed on both the laminated and locally delaminated regions. 
Hence, 


W(a) - WLAM(a) + W LD (a) {4) 

For any region, the work performed may be defined in terms of the 
strain and stress tensors resulting from an applied mechanical load 
on that region as 


W = 



(5) 


Using the assumptions and matrix notation of classical laminated 
plate theory [11], eq (5) becomes 

w = wl X [ (e)^ (a]{^ dz (6) 

k=1 


where the prime in ( e ‘)k represents the transpose of the strain 
vector, and K is the number of plies in the region. For a laminate that 
does not undergo bending, and has equal ply thicknesses throughout 
the laminate, eq (6) becomes 

w - Jf i rfrf (7) 

K k=1 

If the laminate in the region is asymmetric, however, bending and 
coupling contributions to the mechanical strain, and hence, to the 
work done on the region, may be included using the expressions given 
in the appendix. The strains in equations (6) and (7) are the strains 
associated with the applied mechanical load, P, for the region and 
are given by 
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where £x is the axial strain in the region due to the applied 
mechanical load, v X y is the laminate Poisson's ratio for the region, 
and rixy.x is the laminate coefficient of mutual influence of the 
second kind for the region. 

For the laminated region, wlam = w, Ilam = l-a, K = N, and ti_AM - 
Nh, where N is the number of plies in the laminate and h is the ply 
thickness. Therefore, for an applied mechanical load, P, 



P 

WhNE L AM 


(9) 


For the locally delaminated region, wld = w, Ild * a, K = N-n, where n 
is the number of cracked plies, and t|_D = (N-n)h. Therefore 



P 

wh(N-n)E|_D 


( 10 ) 


The Poisson's ratio, v X y, and r| X y iX for each region may be calculated 
from laminated plate theory [11]. Because n X y, x accounts for the 
coupling between the applied axial strain and in-plane shear, it will 
be zero for the laminated region if the original laminate is balanced 
and symmetric, but it may be non-zero for the locally delaminated 
region. 

The stress vector, Mk*, in equations (6) and (7) contains the 
ply stresses associated with the applied axial mechanical load, P, 
for either region and is given by 

rf=[Q]k{e}J* (11) 

where [o]k , is the transformed reduced stiffness matrix for the kth 
ply as defined in ref [11]. Substituting eq (11) into eq (7) yields 
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( 12 ) 


where 


w --2*11 i «r 

* k=1 

u“-l(e)“[Ql < (e| k M 


For the model shown in fig.1, wlam = wld = w, Ilam * I-a, I|_d - a, 
Klam = N, Kld - N-n, tLAM - Nh, and tLD - (N-n)h. Combining these 
expressions with eqs (4) and (12) and noting that each region has 
equal ply thickness, h=t/K, yields 


W(a) = 2wlh u l m am + 2wah [u“ - u" J (14 ) 

where “lam - X u k* ( 15 ) 

k=l 

N-n 

and u LD = X u k* ( 16 ) 

k=1 


Recalling that dA = mwda for a through-width local delamination 
growing uniformly from a matrix ply crack [4,8], and differentiating 
eq (14) with respect to delamination length, a, yields 

' (^) < u “d - u l M am> (17) 

The strain energy released in the laminate as the delamination 
grows is equal to the summation of the strain energy released in the 
laminated and locally delaminated regions. Hence, 

U(a) - U L AM(a) + U LD (a) (18) 

For any region, the strain energy released may be defined in terms of 
the strain and stress tensors for that region as 


8 



dv 


(19) 


Using the assumptions and matrix notation of classical laminated 
plate theory [11], eq (19) becomes 

U -JKl X [ 1 (e)k (Ojk dz (20) 

For a laminate that does not undergo bending, and has equal ply 
thicknesses throughout the laminate, eq (20) becomes 

u -Jttlt £ (e)k(a) k (21) 

2K k.l 

If the laminate in the locally delaminated region is asymmetric, 
however, bending and coupling contributions to the strain, and hence, 
to the strain energy in the region, may be included using the 
expressions given in the appendix. 

The strains in equations (20) and (21) are the total strains 
that contribute to the release of strain energy as the delamination 
grows. This total strain may consist of strains resulting from the 
applied mechanical load (M), the residual thermal stresses in the 
region (T), and the hygroscopic stresses that develop as a function 
of absorbed moisture (H). The total strain may be written as 

{ £ }k - tf + {e}J + (22) 

where the mechanical strain is given by eq (8), with £x* given by eqs 
(9) or (10) for the appropriate region. The thermal strains for any 
region are given by [12] 

( a* - «x \ 

a y -a) AT (23) 

otxy - a k xy | 
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where a- the coefficient of thermal expansion of the region 

a k * the coefficient of thermal expansion of the kth ply in 
the region 

AT « the difference between the stress free temperature 
and the temperature of the ambient environment 

The hygroscopic strains for any region are given by [12] 



| Px-Px 
1 _Py - Py 
, Pxy — Pxy 


AH 


(24) 


where P = the coefficient of moisture expansion of the region 

P k = the coefficient of moisture expansion of the kth ply 
in the region 

AH - the percent weight difference between the ambient 
(wet) condition and the dry condition 

The stress vector, {o}k, in equations (20) and (21) contains the 
ply Stresses associated with the total strain for any region and is 
given by [11] 


(°>k-[a]<(e)k (25) 

Substituting eq (25) into eq (21) and noting that each region has 
equal ply thickness, h=t/K, yields 


K 

U = wlh £ u k 

k=1 


(26) 


where 


u k «i-{e'}k[Q]k {ek 


(27) 
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For the model shown in fig.1, wlam = wld - w, Ilam * l-a, Ild ■ a, 
Klam * N, Kld - N-n, tLAM - Nh, and t|_D - (N-n)h. Combining these 
expressions with eqs (18) and (26) yields 


U(a) = wlh u LAM + wah [u LD - u LAM ] 

N 

where u LAM = X u k 

k=1 


(28) 

(29) 


N-n 

and u LD * X u k 

k=1 


(30) 


Recalling that dA = mwda for a through-width local delamination 
growing uniformly from a matrix ply crack [4,8], and differentiating 
eq (28) with respect to delamination length, a, yields 

^-(£)<Uld-Ulam) (31) 

Substituting eqs (17) and (31) into eq (3) yields the strain energy 
release rate 


G - (£) |2(u l m d - u“ AM ) - (u LD - u LAM )) (32) 

If only the mechanical strain due to the applied load is assumed to 
contribute to the strain energy released in the laminated and 

delaminated regions, then u ld = u ld, u lam = u lam, and 
eq (32) becomes 


G = 




(33) 


Eq (33) will yield the same strain energy release rate as the closed 
form eq (2) derived based on the change in compliance with 
delamination growth. If the laminate is completely dry, only the 
mechanical strain due to the applied load and the strains due to 
residual thermal stresses contribute to the strain energy released in 



the laminated and delaminated regions. For this mechanical plus 
thermal case, the strain energy release rate, G M+T , is given by eq 

|_| 

(32), where AH, and hence ( e )k in eqs (22) and (24), are zero. If the 
laminate is not completely dry, which is typical of polymer matrix 
composites that absorbed moisture from the air in the ambient 
environment [8,12], all of the strains in eq (22) contribute to the 
strain energy released in the laminated and locally delaminated 
regions. For this mechanical plus thermal plus hygroscopic case, the 

strain energy release rate, g m+t+h , is given by eq (32), where AH, 
|_| 

and hence ( e ^k in eqs (22) and (24), are not zero. 


ANALYTICAL RESULTS 

In this section, four laminates will be analyzed to illustrate 
the influence of residual thermal and moisture stresses on G for 
local delamination from matrix cracks. These results simply 
illustrate how eq (32) may be used to quantify G for a combined 
mechanical, thermal, and moisture loading condition. These results 
are not intended as a sensitivity study, nor should they be compared 
to each other, since all three load components are coupled, and only 
their relative magnitudes at the onset of delamination under several 
loadings and environmental conditions can reveal their relative 
significance. 

In the derivation of G M+T+H for edge delamination in ref [12], 
it was noted that the sublaminates in the delaminated region would 
be constrained to the laminate's free thermal and moisture 
expansions in the applied load direction by the test machine grips. 

Therefore, a x and P x for the sublaminates were set equal to a x and P x 
for the original laminate in eqs (23) & (24). The local delaminations 
modeled in this study do not extend to the grip. However, the 
sublaminates in the locally delaminated region may be constrained 
by the sublaminate in the laminated region. Therefore, results for 
particular laminates calculated using eq (32) will be plotted with 
and without this constraint to illustrate the influence of 
constraining the sublaminates' free thermal expansion on the 
calculated strain energy release rate. 

Figure 2 shows the influence of residual thermal and moisture 
stresses on the strain energy release rate associated with a local 
delamination growing from a 90 degree matrix crack in the 0/90 
interfaces of a (02/904) s glass epoxy laminate. Lamina properties 
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from ref [6] were used in the analysis along with the following 
thermal coefficients assumed to be typical for glass epoxy 

ai ■ 3.80 pe/°C (2.11 fie/°F) 

a 2 = 16.7 pe/°C (9.28 pe/°F) 

and the following moisture expansion coefficients that were 
determined for graphite epoxy [12], but were also assumed to be 
reasonable for glass epoxy 

pi = 0 pe/weight% 

p 2 = 5560 pe/weight% 

A mechanical axial load per unit width, N x , of 648 kN/m (3703 
Ibs/in.), corresponding to an axial strain of 0.01, and a AT of -100°C 
(-180°F) were applied. This AT corresponds to the difference 
between the room temperature condition of 21 °C (70°F) and a 121°C 
(250° F) cure temperature for glass epoxy. The glass epoxy ply 
thickness was assumed to be 0.203 mm (0.008 in.). Strain energy 
release rates are plotted as a function of the percentage moisture 
weight gain assuming either a constrained or unconstrained 
sublaminate as discussed previously. The strain energy release rate 
caused by mechanical loading only, G M , calculated from eq (33) 
agreed with the value calculated from eq (2) and is shown on the 
ordinate. The strain energy release rate corresponding to a 
mechanical and thermal loading only, G M+T , is also shown on the 
ordinate. 

For the constrained sublaminate case, G M+T , as calculated by 
eq (32), is higher than G M for the same applied load. However, if the 
laminate absorbs moisture, the residual thermal stresses are 
relaxed, and the strain energy release rate, G M+T+H , decreases 
depending on the percentage of moisture weight gain, AH. 

For the constrained sublaminate case shown in fig. 2, a moisture 
weight gain of approximately 0.2% completely relaxes the residual 
thermal stresses, resulting in a strain energy release rate, G M+T+H , 
that is equal to G^ due to the mechanical load alone. For the 
unconstrained sublaminate case, however, G M+T calculated by eq 
(32) is only slightly higher than G M for the same applied load. 
Furthermore, as the laminate absorbs moisture, the strain energy 
release rate, G M+T+H , decreases to G M at AH ■ 0.2%, but then 
increases monotonically with increasing AH. 


13 


Fig. 3 shows a similar relationship for the strain energy 
release rate associated with a local delamination growing from a 
90 degree matrix crack in the -45/90 interfaces of a (45/-45/904) s 
glass epoxy laminate subjected to an N x of 414 kN/m (2363 Ibs/in.), 
corresponding to an axial strain of 0.01, and a AT of -100°C (-180°F). 
For the constrained sublaminate case, the strain energy release rate 
caused by mechanical loading only, G M , increases slightly when the 
residual thermal stresses are included, G M+T , and decreases when 
the moisture contribution is included, G M+T+H . For the unconstrained 
sublaminate case, however, G M+T is nearly identical to G M for the 
same applied load. Furthermore, as the laminate absorbs moisture, 
the strain energy release rate, G M+T+H , increases very slightly with 
increasing AH. 

In ref. [8], a strain energy release rate analysis was 
developed for local delaminations forming in the 0/-0 interfaces of 
a (O/0/-0 ) s laminate from a -0 degree matrix crack, where 0 is 
greater than zero but less than ninety degrees. For these cases, the 
sublaminates that comprise the locally delaminated (LD) region 
consist of a pair of asymmetric (0/0) sublaminates. Under the 
applied tension load, these sublaminates may undergo both shear and 
bending deformations depending on the degree of shear-extension 
and bending-extension coupling present. To account for these 
coupling effects, the modulus of the locally delaminated region, Eld, 
in eq (2) was determined by prescribing an applied mechanical 
strain, e x , setting the applied N y and M y equal to 0, and setting y xy , k x 
and K xy equal to 0. The N x that results from this analysis was used 
to determine Eld as 


ELD - N x /tLDE X (34) 

which may be used in eq (2) to calculate G M and in eq (10) to 

calculate (eiCW Furthermore, This loading may be prescribed to 
determine the ply stresses and strains used to calculate the strain 

M M 

energy and work terms in eq (32) where u lam, u ld, u lam, and u ld are 
described in the appendix. 

Figure 4 shows the influence of residual thermal and moisture 
stresses on the strain energy release rate associated with a local 
delamination growing from a -15 degree matrix crack in the 15/-15 
interfaces of a (0/15/-15)$ graphite epoxy laminate. Lamina 
properties from ref [8] were used in the analysis along with the 
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following thermal and moisture coefficients for graphite epoxy from 
ref [12] 


<xi = -0.41 pe/°C (-0.23 pe/°F) 

02 = 26.8 pe/°C (14.9 pe/°F) 

Pi = 0 pe/weight% 

P 2 = 5560 pe/weight% 

A mechanical axial load per unit width, N x , of 1400 kN/m (8000 
Ibs/in.) and a AT of -156°C (-280°F) were applied. This AT 
corresponded to the difference between the room temperature 
condition of 21 °C (70°F) and a 177°C (350°F) cure temperature for 
graphite epoxy. The N x chosen is close to the average N x at the onset 
of damage in (O 2 / 1 52/-1 52) s AS4/3501-6 graphite epoxy laminates 
[7,8]. In the analysis, the graphite epoxy ply thickness was assumed 

to be 0.254 mm (0.01 in.) which is twice the measured ply thickness. 

This allowed the analysis of a 6-ply (0/15/-15) s laminate instead 
of a 12-ply (02/152/-152 )s laminate. Strain energy release rates are 
plotted as a function of the percentage moisture weight gain 
assuming either an unconstrained sublaminate or a sublaminate 
whose free thermal and moisture expansions in the load (X) direction 
were constrained, as discussed previously. Furthermore, because 
these sublaminates experience shear-extension coupling, G was also 
calculated for sublaminates whose free thermal and moisture shear 

expansions were constrained. Hence, a xy and P xy for the 

sublaminates were set equal to a xy and P xy for the original laminate. 
G was calculated with either a shear constraint alone (XY) or a shear 
constraint in addition to the constraint in the load direction (X&XY). 

In all cases, the strain energy release rate caused by 

mechanical loading only, G M , calculated from eq (33) agreed with the 
value calculated from eq (2) and is shown on the ordinate. The strain 
energy release rate corresponding to a mechanical and thermal 
loading only, G M + T , is also shown on the ordinate. For the 
unconstrained and XY constrained cases, G M+T calculated from eq 
(32), is higher than G M for the same applied load. However, for the X 
constrained and X&XY constrained cases, G M+T calculated from eq 
(32) is nearly identical to G M for the same applied load. For the 
unconstrained and XY constrained cases, as the laminate absorbs 
moisture, the residual thermal stresses are relaxed, and the strain 
energy release rate, G M+T+H calculated from eq (32) decreases 
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with increasing moisture weight gain, AH. However, for the X 
constrained and X&XY constrained cases, G M+T+H calculated from eq 
(32) deviates only slightly from G M with increasing moisture weight 
gain, AH. 

Fig. 5 shows a similar relationship for the strain energy 
release rate associated with a local delamination growing from a 
-45 degree matrix crack in the 45/-45 interfaces of a (0/45/-45) s 
graphite epoxy laminate. An N x of 1400 kN/m (8000 Ibs/in.) and a AT 
of -156°C (-280° F) were applied. In all cases, the strain energy 
release rate caused by mechanical loading only, G M , calculated from 
eq (33) agreed with the value calculated from eq (2) and is shown on 
the ordinate. The strain energy release rate corresponding to a 
mechanical and thermal loading only, G^ + T, is also shown on the 
ordinate. For the unconstrained and XY constrained cases, G M + T 
calculated from eq (32) is significantly higher than G^ for the same 
applied load. However, for the X constrained and X&XY constrained 
cases, G M+ T calculated from eq (32) is only slightly higher than G M 
for the same applied load. For the unconstrained and XY constrained 
cases, as the laminate absorbs moisture, the residual thermal 
stresses are relaxed, and the strain energy release rate, g m+t+ h 
calculated from eq (32) decreases with increasing moisture weight 
gain, AH. However, for the X constrained and X&XY constrained cases, 
qM+t+h calculated from eq (32) deviates only slightly from G M with 
increasing moisture weight gain, AH. 

DISCUSSION 

For the local delaminations that originate from 90 degree 
matrix ply cracks (figures 2,3), the constrained values for GM+t+h 
are probably more physically realistic, because these delaminations 
are localized along the laminate's length. Hence, it is possible to 
visualize how the undamaged sublaminate may constrain the free 
thermal expansion of the locally damaged region in the load 
direction. These constrained results for the local delamination from 
90 degree matrix cracks are also consistent with the constrained 
results for edge delamination [12] where G M+T is greater than G M , 
but G m+t+h decreases with increasing moisture weight gain. 

For the local delaminations that originate from angle ply (non- 
90 degree) matrix ply cracks (figures 4,5), the XY constraint has 
very little effect on g m+t+h when the free thermal and moisture 
expansions in the X direction are either constrained or 
unconstrained. However, the X constrained values of G M+T+H are 
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significantly different than the unconstrained values. Hence, the X or 
X&XY constrained values for G M+T+H are probably more physically 
realistic because these delaminations are also localized along the 
laminate's length and extend only slightly into the laminate's width 
[7,8]. However, these constrained results indicate that residual 
thermal and moisture contributions to G for local delaminations that 
originate from -0 degree ply cracks in (O/0/-0) s graphite epoxy 
laminates are very small, unlike the more significant effects 
observed for edge delamination in quasi-isotropic graphite epoxy 
laminates [12]. 

Tests on laminates that are dry could conceivably be compared 
to tests that have absorbed a known moisture content to quantify 
the significance of these thermal and moisture contributions. This 
experimental verification has been attempted for edge delamination 
with some success [12]. However, as noted in ref [7], under 
monotonic loading of (0/1 5/-1 5) s graphite epoxy laminates, the 
first indication of damage that can be documented shows both 
matrix cracking and local delaminations present. Hence, the observed 
dependence of moisture content on the load or strain at the onset of 
damage may reflect the influence of thermal and moisture stresses 
on matrix cracking more than their influence on local delamination. 
Similar tests were conducted under cyclic loading [8], where it was 
documented that matrix cracking preceded the onset of local 
delamination. Although the moisture content before testing was 
determined in ref.8, no attempt was made to test dry laminates in 
fatigue because of their tendency to absorb moisture from the air 
during the course of the fatigue test. 

Clearly, a precise experimental verification of eq (32) is 
difficult to achieve. However, the accuracy of eq (32) could be 
verified by comparing to 3D finite element analysis results that also 
incorporate thermal and moisture effects. Hopefully, these 
comparisons should help confirm the correct constraint assumptions 
for the cases studied. 

CONCLUSIONS 

A laminated plate theory analysis was developed to calculate 
the strain energy release rate associated with local delaminations 
originating at matrix ply cracks in laminates subjected to axial 
loads. The analysis included the contribution of residual thermal and 
moisture stresses to the strain energy released. Examples were 
calculated for the strain energy release rate associated with local 
delaminations originating at 90 degree and angle ply (non-90 degree) 
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matrix ply cracks in glass epoxy and graphite epoxy laminates, 
respectively. The solution developed may be used to assess the 
relative contribution of mechanical, residual thermal, and moisture 
stresses on the strain energy release rate for local delamination for 
a variety of layups and loading combinations. 
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APPENDIX - COUPLING EFFECTS 


In appendix A of ref. 12, expressions were derived for the 
strain energy associated with an asymmetric sublaminate that 
experienced bending-extension coupling. For these asymmetric 
sublaminates, the total strain at any point in the sublaminate is 
given by 


(e) = je°| + z{k} (A1 ) 

where ( E °l are the midplane strains, are the midplane curvatures, 
and z is the distance from the midplane to the point of interest. The 
strain energy for these sublaminates was derived in ref. 12 as 

K 

U = wl X ,k Uk 

k=i (A2) 

where u k “ u ek + u ck + u bk (A3) 

The individual terms in eq (A3) correspond to the strain energy in 
the asymmetric sublaminate corresponding to extension, u ek, 
bending, u bk, and bending-extension coupling, u ck, as defined in eqs 
(A10-12) in ref. 12. Unfortunately, there were several typographical 
errors in these equations. On page 45 of appendix A in ref. 12, in the 
introduction to eq (A8), the "mechanical" strain is defined as the 
"total" strain less the free thermal and moisture strains of the ply. 
This "mechanical" strain is 

( E °| - (a) k AT - (p) k AH 

where 

{e 0 } = ( C °) M + {a}AT + |p)AH 

Hence this "mechanical" strain was not the mechanical term i e ”l 
alone, but wasTactually the total strain described in eqs (22-24) of 
this paper, and eqs (12-14) of ref. 12, which includes the mechanical 
strain and a thermal and moisture contribution corresponding to the 
difference in the free thermal and moisture expansions or 


(A4) 


(A5) 

l nlM 
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contractions of the laminated regions minus the thermal and 
moisture expansions or contractions of the individual plies in the 
regions. These are the strains that contribute to residual thermal 
and moisture stresses, and hence, contribute to the strain energy 
that is stored in each laminated region. In eq (A10) of appendix A in 
ref.12, This strain should^ have appeared on both sides, not just on 

the right hand side, of [ok Hence, eq (A10) in ref.12 should have 
been 


- («) k AT - (p) k AHHa]k({c°| - (cc)kAT - (p) k AH) 


(A6) 


Equation (All) was shown correctly in ref.12 as 

Ubk=i(z^ + ZkZk-1 + z?k i){k}'[qL {*} 

6 


(A7) 


However, eq (A12) in ref.12 should have been 

Uck = i{Zk + Zk i){k}'[q ] k ( ft 0 ) - {a} k AT - (p) k AH) 


(A8) 


The calculations and results reported in ref.12 were performed 
correctly using eqs (A6-A8) above, even though the errors noted 
were made in printing the text. 

For the asymmetric locally delaminated regions modeled in 

,.M 

this paper, u k in eqs (15 &16) and Uk in eqs (29 & 30) were 
calculated using equation A3 (where u ek- u bk and u ck are given by 
equations A6,A7, &A8) instead of using equations 13 and 27 alone. 

n M 

These u k and Uk values were used in eqs (15,16,29,&30) to determine 
M u 

u lam, u ld, u lam, and u ld, respectively. 
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Figure 1. Model of local delamination. 


Fig. 2 q M+T+H FOR LOCAL DELAMINATION 
FROM 90 DEGREE MATRIX CRACK 
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Fig. 3 q M+T+H FOR LOCAL DELAMINATION 
FROM 90 DEGREE MATRIX CRACK 



25 






Fig. 4 g M+T+H FOR LOCAL DELAMINATION 
FROM -15 DEGREE MATRIX CRACK 
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Fig. 5 g M+T+H FOR LOCAL DELAMINATION 
FROM -45 DEGREE MATRIX CRACK 
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